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ABSTRACT 

Aims. We seek to identify the primaiy agents causing Forbush decreases (FDs) observed at the Earth in high rigidity cosmic rays. In 

particular, we ask if such FDs are caused mainly by coronal mass ejections (CMEs) from the Sun that are directed towards the Earth, 

or by their associated shocks. 

Methods. We use the muon data at cutoff rigidities ranging from 14 to 24 GV from the GRAPES-3 tracking muon telescope to 

identify FD events. We select those ED events that have a reasonably clean profile, and can be reasonably well associated with an 

Earth-directed CME and its associated shock. We employ two models: one that considers the CME as the sole cause of the FD (the 

CME-only model) and one that considers the shock as the only agent causing the FD (the shock-only model). We use an extensive 

set of observationally determined parameters for both these models. The only free parameter in these models is the level of MHD 

turbulence in the sheath region, which mediates cosmic ray diffusion (into the CME, for the CME-only model and across the shock 

sheath, for the shock-only model). 

Results. We find that good fits to the GRAPES-3 multi-rigidity data using the CME-only model require turbulence levels in the CME 

sheath region that are only slightly higher than those estimated for the quiet solar wind. On the other hand, reasonable model fits with 

the shock-only model require turbulence levels in the sheath region that are an order of magnitude higher than those in the quiet solar 

wind. 

Conclusions. This observation naturally leads to the conclusion that the Earth-directed CMEs are the primary contributors to FDs 

observed in high rigidity cosmic rays. 

Key words. Forbush decrease. Coronal mass ejection(CME), Cosmic rays 



> 

m 

o 
m 



X 



1. Introduction 

Forbush decreases are short-term decreases in the intensity of 
the galactic cosmic rays at the Earth. They are typically caused 
by the effects of interplanetary counterparts of coronal mass 
ejections (CMEs) from the Sun, and also the corotating inter- 
action regions (CIRs) between the fast and slow solar wind 
streams from the Sun. We concentrate on CME driven Forbush 
decreases in this paper. The near-Earth manifestations of CMEs 
from the Sun typically have two major components: the inter- 
planetary counterpart of the CME, commonly called an ICME, 
and the shock which is driven ahead of it. ICMEs which pos- 
sess certain well-defined criteria such as plasma temperature 
depressions and smooth magnetic field rotations (e.g., Burlaga 
et al. 1981, Bothmer & Schwenn 1998) are called magnetic 
clouds, while others are often classified as "ejecta". The rela- 
tive contributions of shocks and ICMEs in causing Forbush de- 
creases is a matter of debate. For instance, Zhang & Burlaga 
(1988), Lockwood, Webber, & Debrunner (1991) and Reames, 
Kahler & TyUca (2009) argue against the contribution of mag- 
netic clouds to Forbush decreases. On the other hand, other 
studies (e.g., Badruddin, Yadav, & Yadav, 1986; Sanderson et 
al.,1990; Kuwabara et al 2009) concluded that magnetic clouds 
can make an important contribution to FDs. There have been 
recent conclusive associations of Forbush decreases with Earth- 
directed CMEs (Blanco et al 2013; Oh & Yi 2012). Cane (2000) 



introduced the concept of a "2-step" FD, where the first step 
of the decrease is due to the shock and the second one is due 
to the ICME. Based on an extensive study of ICME-associated 
Forbush decreases at cosmic ray energies between 0.5 - 450 
MeV, Richardson and Cane (2011) conclude that shock and 
ICME effects are equally responsible for the Forbush decrease. 
They also find that ICMEs that can be classified as magnetic 
clouds are usually involved in the largest of the Forbush de- 
creases they studied. From now on, we will use the term "CME" 
to denote the CME near the Sun, as well as its counterpart ob- 
served at the Earth. 

In this paper we have used Forbush decrease data from 
GRAPES-3 tracking muon telescope located at Ooty (11.4°N 
latitude, 76.7°E longitude, and 2200 m altitude) in south India. 
The GRAPES-3 muon telescope records the flux of muons in 
nine independent directions (labeled NW, N, NE, W, V, E, SW, 
S and SE), and the geomagnetic cutoff" rigidity over this field 
of view varies from 12 to 42 GV. The details of this telescope 
are discussed in Hayashi et al. (2005), Nonaka et al. (2006) and 
Subramanian et al (2009). Thus, the GRAPES-3 telescope ob- 
serves the cosmic ray muon flux in nine independent directions 
with varying cut-off" rigidities simultaneously. The high muon 
counting rate measured by the GRAPES-3 telescope results in 
extremely small statistical errors, allowing small changes in the 
intensity of the cosmic ray flux to be measured with high preci- 
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CME Sheath 



Fig. 1. A schematic of the CME-shock system. The CME is 
modeled as a flux rope structure. The undulating lines ahead of 
the shock denote MHD turbulence driven by the shock, while 
those in the CME sheath region denote turbulence in that region. 



2. Event shortlisting criteria 

2. 1 . First and second shortlists: cliaracteristics of tlie 
Forbusli decrease 

We have examined all Forbush decrease events observed by the 
GRAPES-3 muon telescope from 2001 to 2004. We then short- 
listed events that have magnitudes > 0.25 % and possess a 'rela- 
tively clean' profile comprising a sudden decrease followed by a 
gradual exponential recovery. While the figure of 0.25 % might 
seem rather small by neutron monitor standards, we emphasize 
that the largest events obsei-ved with GRAPES-3 have magni- 
tudes of ~ 1 %. This list, which contains 72 events, is called 
shortlist 1 . We next correlate the events in shortlist 1 with lists 
of magnetic clouds near the Earth observed by the WIND and 
ACE spacecrafts (Huttunen et al 2005; Lynch et al. 2003; A. 
Lara, private communication) and select only those that can be 
reasonably well connected with a near-Earth magnetic cloud and 
this set is labeled shortlist 2 (Table [T]i- The decrease minimum 
for most of the FD events in shortlist 2 lie between the start and 
the end of the near-Earth magnetic cloud. 



sion. Thus a small drop (~0.2%) in the cosmic ray flux during 
a Forbush decrease event can be reliably detected. This is possi- 
ble even in the presence of the diurnal anisotropy of much larger 
magnitude (~L0%), through a filtering technique described in 
Subramanian et al (2009) (referred to hereafter as S09). 

A schematic of the CME, which is assumed to have a flux- 
rope geometry (Vourlidas et al 2012) together with the shock it 
drives, is shown in Figure [1] The shock drives turbulence ahead 
of it, and there is also turbulence in the CME sheath region (e.g., 
Manoharan et al 2000; Richardson & Cane 2011). 

Instead of treating the entire system shown in Figure [1] 
which would be rather involved, we consider two separate mod- 
els. The first, that we call the "CME-only" model, is one where 
the Forbush decrease is assumed to be exclusively due to the 
CME, which is progressively populated by high energy cosmic 
rays as it propagates from the Sun to the Earth. The preliminary 
idea behind this model was first sketched by Cane, Richardson 
& Wibberenz (1995) and developed in detail in S09. The work 
described here addresses multi-rigidity data, which is a major 
improvement over S09. We will describe several other salient 
improvements in the CME-only model in subsequent sections. 
The second model, which we call the "shock-only" model, is 
one where the Forbush decrease is assumed arise only due to a 
propagating diffusive barrier, which is the shock driven by the 
CME (e.g., Wibberenz et al 1998). The diffusive barrier would 
act as a shield for the galactic cosmic ray flux, resulting in a 
lower cosmic ray density behind it. In treating these two models 
separately, we aim to identify which is the dominant contributor 
to the observed Forbush decrease; the CME, or the shock. 

We identify Forbush decreases in the GRAPES-3 data that 
are associated with both near-Earth magnetic clouds and possi- 
bly shocks driven by them. We describe our event shortlisting 
criteria in the next section. We then test the extent to which each 
of the models (the CME-only and the shock-only model) satis- 
fies the multi-rigidity Forbush decrease data from the GRAPES- 
3 muon telescope. In the subsequent analysis the use of cutoff 
rather than median rigidity is preferred for the following reason. 
The cutoff rigidity in a given direction represents the threshold 
rigidity of incoming cosmic rays and the magnitude of Forbush 
decrease is a sensitive function of it. On the other hand the me- 
dian rigidity is comparatively insensitive to the magnitude of the 
Forbush decrease. 



2.2. Third shortlist: CME velocity profile 

Since we are looking for Forbush decrease events that are associ- 
ated with shocks as well as CMEs, we examine the CME catalog 
[littp : / /cdaw.gsfc.nasa.gov/CMEJist/) for a near-Sun CME 
that can reasonably correspond to the near-Earth magnetic cloud 
that we associated with the Forbush decrease in forming shortlist 
2. In S09 it was assumed that the CME propagates with a con- 
stant speed from the Sun to the Earth. In this paper, we adopt a 
more realistic, two-step velocity profile that is described below. 
The data from the LASCO coronograph aboard the SOHO 
spacecraft {http : / /cdaw.gsfc.nasa.gov/CMEJist/) provide 
details about CME propagation up to a distance of ^ 30 Rq from 
the Sun. We fit the following velocity profile to the LASCO data 
points: 



Vi + ait , for R(t) < Rn 



(1) 



where Vi and ai are the initial velocity and acceleration of CME 
respectively, and R,,, is the heliocentric distance at which the 
CME is last observed in the LASCO field of view. For distances 

> Rm, we assume that the CME dynamics are governed exclu- 
sively by the aerodynamic drag it experiences due to momentum 
coupling with the ambient solar wind. For heliocentric distances 

> Rm, we therefore use the following widely used ID differential 
equation to determine the CME velocity profile: (e.g., Borgazzi 
et al 2009) 



mcMEVi— — - -CdPswAcme(V2 - Vs„) , 
oR 2 



R(t) > R„ 



(2) 



where mcME is the CME mass, Co is the dimensionless drag 
coeflicient, psw is the solar wind density. Acme = '^'^cme ^^ ^^^ 
cross-sectional area of the CME and Vsw is the solar wind speed. 
The boundary condition used is V2 = Vn, at R(t) - Rm. The CME 
mass mcME is assumed to be lO'^ g and the solar wind speed 
Vsw is taken to be equal to 450 km/s. The solar wind density psw 
is given by the model of Leblanc et al (1998). The composite 
velocity profile for the CME is defined by 



VcME 






R(t) < R„ 
R(t) > R„ 



(3) 



The total travel time for the CME is J dR/VcME, where Ri is 
the heliocentric radius at which the CME is first detected and Rf 
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Table 1. Selected Events 



Event 


Magnitude(Ver) (%) 


FD onset (UT) 


FD min (UT) 


MC start (UT) 


MC stop (UT) 


CME launch (UT) 


11/04/01 


2.86 


11/04/01, 12:00 


12/04/01, 18:00 


11/04/01,23:00 


12/04/01, 18:00 


10/04/01, 05:30 


17/08/01 


1.03 


16/08/01, 22:34 


18/08/01, 05:00 


18/08/01, 00:00 


18/08/01,21:30 


15/08/0 1, 23:54 


29/09/01 


2.44 


29/09/01, 12:43 


31/09/01, 01:00 


29/09/01, 16:30 


30/09/01, 11:30 


. . . not clear. . . 


24/11/01 


1.56 


24/11/01,03:21 


25/11/01, 15:00 


24/11/01, 17:00 


25/11/01, 13:30 


22/11/01,22:48 


23/05/02 


0.93 


23/05/02,02:10 


23/05/02, 23:00 


23/05/02,21:30 


25/05/02, 18:00 


. . . not clear. . . 


7/09/02 


0.97 


7/09/02, 14:52 


8/09/02, 13:00 


7/09/02, 17:00 


8/09/02, 16:30 


05/09/01, 16:54 


30/09/02 


0.97 


30/09/02, 11:30 


31/09/02, 08:00 


30/09/02, 22:00 


1/09/02, 16:30 


. . . not clear. . . 


20/11/03 


1.16 


20/11/03, 10:48 


24/11/03,04:00 


21/11/03,06:10 


22/11/03,06:50 


18/11/03,08:50 


26/07/04 


2.13 


26/07/04, 15:36 


27/07/04, 11:00 


27/07/04, 02:00 


28/07/04, 00:00 


25/07/04, 14:54 



Notes. Events that can clearly be associated with near-Sun CMEs and agree with the composite velocity profile (§ 2.2) have a CME launch time 

(near the Sun) mentioned against them. These events form our final shortlist (shortlist 3). 

Magnitude(Ver) : Magnitude of Forbush decrease in vertical direction 

FD onset : Time of FD onset in UT 

FD minimum : Time of FD minimum in UT 

MC start : Magnetic cloud start time in UT 

MCstop : Magnetic cloud stop time in UT 

CME near Sun : Time at which CME was first observed in the LASCO FOV 
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Fig. 2. A plot of the velocity profile for the CME correspond- 
ing to the 24 November 2001 FD event. The CME was first ob- 
served in the LASCO FOV on 22 November 2001 . The solid line 
shows the first stage governed by LASCO observations (Eq [T}, 
where the CME is assumed to have a constant deceleration. The 
dashed line shows the second stage (Eq|2|, where the the CME 
is assumed to experience an aerodynamic drag characterized by 
a constant Cd. 



is equal to 1 AU. We have used a constant drag coefficient Cd 
and adjusted its value so that the total travel time thus calculated 
matches the time elapsed between the first detection of the CME 
in the LASCO FOV and its detection as a magnetic cloud near 
the Earth. We have retained only those events for which it is 
possible to find a constant Cd and this criterion is satisfied. We 
have also eliminated magnetic clouds that could be associated 
with multiple CMEs. This defines our final shortlist, which we 
call shortlist 3. The events that have an entry in the last column 
labeled 'CME launch' in Table [1] comprise shoitlist 3. We note 
that we have adjusted the parameter Cd so that the final CME 
speed obtained from Eq (O is close to the observed magnetic 
cloud speed near the Earth. It is usually not possible to find a 
Cd that will yield an exact match for the velocities as well as 
the total travel times (e.g., Lara et al 2011). Figure |2] shows an 
example of the composite velocity profile (given by Eqs[T]and|2]i 
for one such representative CME: the one that was first observed 
in the LASCO FOV on 22 November 2001, and resulted in a FD 
on 24 November 2001. 



3. Models for Forbush decreases 

We apply two different models to the FD events in Table [1] the 
first one is the CME-only model, which assumes that the Forbush 
decrease owes its origin only to the CME. The second one is the 
shock-only model, which assumes that the Forbush decrease is 
exclusively due to the shock, which is approximated as a diffu- 
sive barrier. Although both the shock and the CME are expected 
to contribute to the Forbush decrease, our treatment seeks to de- 
termine which one of them is the dominant contributor at rigidi- 
ties ranging from 14 to 24 GV Before describing the models, we 
discuss the perpendicular diffusion coefficient. 

3. 1 . Perpendicular diffusion coefficient 

We use an isotropic perpendicular diffusion coefficient (Dj^) to 
characterize the penetration of cosmic rays across large scale, or- 
dered magnetic fields. We envisage a CME, which has a flux rope 
structure, propagates outwards from the Sun, driving a shock 
ahead of it (see, e.g., Vourlidas et al 2012). The flux rope CME- 
shock geometry is illustrated in Figure [T] The CME sheath re- 
gion between the CME and the shock is known to be turbu- 
lent (e.g., Manoharan et al 2000) and it is well accepted that 
it has a significant role to play in determining Forbush decreases 
(Badruddin 2002; Yu et al 2010). The isotropic perpendicular 
diffusion coefficient Dj^ characterizes the penetration of cosmic 
rays through the ordered, compressed large-scale magnetic field 
near the shock, as well as across the ordered magnetic field of 
the flux rope CME. In diffusing across the shock, the cosmic rays 
are affected by the turbulence ahead of the shock, and in diffus- 
ing across the magnetic fields bounding the flux rope CME, the 
cosmic rays are affected by the turbulence in the CME sheath 
region. 

The subject of charged particle diffusion across field lines in 
the presence of turbulence is a subject of considerable research. 
Analytical treatments include the so-called "classical" scattering 
theory (e.g., Giacalone & Jokipii (1999) and references therein), 
and the nonlinear guiding center theory (Matthaeus et al 2003; 
Shalchi 2010) for perpendicular diffusion. Numerical treatments 
include Giacalone & Jokipii (1999), Casse et al (2002), Candia 
& Roulet (2004) and Tautz & Shalchi (201 1). For our purposes, 
we seek a concrete, usable prescription for the D^ that can in- 
corporate as many observationally determined quantities as pos- 
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sible. We find that the analytical fits to extensive numerical sim- 
ulations provided by Candia & Roulet (2004) best suit our re- 
quirements. We note, in particular, that their results not only 
reproduce the standard results of Giacalone & Jokipii (1999) 
and Casse et al (2002) but also extend the regime of validity 
to include strong turbulence and high rigidities. Our approach 
is similar to that of Effenberger et al (2012), who adopt em- 
pirical expressions for the perpendicular diffusion coefficients. 
It may be mentioned, however, that they allow for the possi- 
bihty of anisotropic perpendicular diffusion. We treat the case 
of radial diffusion across the largely azimuthal magnetic fields 
bounding a flux rope CME structure, and we therefore need only 
one (isotropic) perpendicular diffusion coefficient (D^). 

In the formulation of Candia & Roulet (2004), the perpen- 
dicular diffusion coefficient D^ is a function of the quantity p 
(which is closely related to the rigidity and indicates how tightly 
the proton is bound to the magnetic field) and the level of turbu- 
lence o"^. Our characterization of T)±{p, cr^) follows that of S09, 
who adopt the analytical fits to Monte Carlo simulations of par- 
ticle transport in turbulent magnetic fields given by Candia & 
Roulet (2004). The quantity p is related to the rigidity Rg by 

where Ri is the Larmor radius, Bq is the strength of the rel- 
evant large-scale magnetic field. For the CME-only model, Bq 
refers to the large-scale magnetic field bounding the CME, and 
for the shock-only model, it refers to the enhanced large-scale 
magnetic field at the shock. In writing second step in Eq (HJi, we 
have related the Larmor radius to the rigidity Rg by 



R,,(.) = SS, 

£>0 



(5) 



For the CME-only model, we adopt L^ax = 2 R(T), where 
R(T) is the radius of the near-Earth magnetic cloud. This is in 
contrast with S09, where L,nax was taken to be lO*" km, which is 
the approximate value for the outer scale of the turbulent cascade 
in the solar wind. For the shock-only model, on the other hand, 
we assume that L,nax is equal to 1 AU. 

The turbulence level cr^ is defined (as in Candia & Roulet 
[2004] and S09) to be 



CT^^ 



(Br') 

Bo' 



(6) 



where Bi is the fluctuating part of the turbulent magnetic 
field and the angular braces denote an ensemble average. 

For the sake of completeness, we reproduce the full expres- 
sion for the isotropic perpendicular diffusion coefficient we use 
in this work. It is the same as the one used in S09, and is taken 
from Candia & Roulet (2004). 

The diffusion coefficient due to scattering of particles along 
the mean magnetic field D\\ is given by 



Dn^cL 



max P 




2(1 -y) 



(7) 



where c is the speed of light and the quantities N\\, y and py are 
constants specific for different kinds of turbulence. The perpen- 
dicular diffusion coefficient (D^) is related to the parallel one 

(D||) by 



£>ii 



[NAcr'^r^, 



^Ao-y^\iz 



p<0.2 
p>0.2 



(8) 



The quantities A^^ and aj_ are constants specific to differ- 
ent kinds of turbulent spectra. In this work, we assume the 
Kolmogorov turbulence spectrum in our calculations. We use A^y 
= 1.7 , y = 5/3 , pii =0.20 , A^^ =0.025 and a^= 1.36 ( Table 1, 
Candia & Roulet [2004] ). Eq Q together with Eq © defines 
the isotropic perpendicular diffusion coefficient we use in this 
work. 



3.2. CME-only model 

The basic features of the CME-only model are similar to that 
used in S09 and here we only highlight the areas where there are 
significant differences from S09. There are practically no high- 
energy galactic cosmic rays inside the CME when it starts out 
near the Sun. The cosmic rays diffuse into it from the surround- 
ings via perpendicular diffusion across the closed magnetic field 
lines as it propagates through the heliosphere as shown schemat- 
ically in Figure |3] Near the Earth, the difference between the 
(relatively lower) cosmic ray proton density inside the CME and 
that in the ambient medium appears as the Forbush decrease. We 
proceed to obtain an estimate of the cosmic ray proton density 
inside the CME produced by the cumulative effect of diffusion. 
The flux F of protons diffusing into the CME at a given time 
depends on the isotropic perpendicular diffusion coefficient Dj^ 
and the density gradient 5Na/cJr, and can be written as 



F (cm 



OT 



(9) 



As mentioned earlier, the D^ characterizes diffusion across the 
(largely closed) magnetic fields bounding the CME and Na is 
the ambient density of high energy protons. The total number of 
cosmic ray protons that will have diffused into the CME after a 
time T is related to the diffusing flux by 



Ui 



XT fT ^xT 

A(t)F(t)dt=J D±A(t)^ 



dt 



(10) 



where A(t) is the cross-sectional area of the CME at a given 
time t. According to our convention, the CME is first observed in 
the LASCO FOV at t = and it reaches the Earth at t = T. The 
ambient density gradient cJNa/i9r is approximated by the follow- 
ing expression, which is significantly different from that used in 
S09: 



5Na Na 



(11) 



where L is the gradient lengthscale. Observations of the density 
gradient lengthscale L exist only for a few rigidities. We use the 
observations of Heber et al (2008), who quote a value of L"' = 
4.7 % AU"' for 1.2 GV protons. We take this as our reference 
value. In order to calculate L for other rigidities (in the 14-24 
GV range that we use here), we assume that L oc R^' . This is 
broadly consistent with the observation (De Simone et al, 201 1) 
that the density gradient lengthscale is only weakly dependent on 
rigidity. For a given rigidity, we also need the value of L all the 
way from the Sun to the Earth. In order to do this, we recognize 
that L near the CME/magnetic cloud will not be the same as its 
value in the ambient solar wind. We adopt L (x BMc/Ba, where 
Bmc is the large-scale magnetic field bounding the CME and B-^ 
is the (weaker) magnetic field in the ambient medium outside 
the CME. Furthermore, while Bmc varies according to Eq ([T2l i 
below, the ambient field B^ of the Parker spiral in the ecliptic 
plane varies inversely with heliocentric distance. 
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Cosmic rays 




Fig. 3. A cartoon illustrating a flux rope CME expanding and propagating away from the Sun. High energy galactic cosmic rays 
diffuse into the CME across its bounding magnetic field. 



We assume that the magnetic flux associated with the CME is 
"frozen-in" with it as it propagates. In other words, the product 
of the CME magnetic field and the CME cross-sectional area 
remains constant. This is a fairly well-founded assumption (e.g., 
Kumar & Rust 1996; Subramanian & Vourlidas 2007). One can 
therefore relate the CME magnetic field B()(t) at a given time t, 
to the value Bmc measured in the near-Earth magnetic cloud by 



Bo(t) - Bmc 



R(T) 



R(t) 



(12) 



where R(T) is the radius of the magnetic cloud observed at the 
Earth and R{t) is its radius at any other time t during its passage 
from the Sun to the Earth. The CME radius R(t) and R(T) are 
related via Eq (fTsT i below. We emphasize that the magnetic field 
referred by Eq (ITZt refers to the magnetic field bounding the 
CME, and not the ambient magnetic field outside it. 

We model the CME as an expanding cylindrical flux rope 
whose length increases with time as it propagates outwards. Its 
cross-sectional area at time t is 



A(t) = 27rL(t)R(t) 



(13) 



where L(t) is the length of the flux-rope cylinder at time t, and is 
related to the height H(t) of the CME above the solar limb via 



L(t) = ;rH(t). 



(14) 



We note that Eq (fT4l) differs from the definition used in 
S09 by a factor of 2. We assume that the CMEs expand in 
a self-similar manner as they propagate outwards. 3D flux 



rope fittings to CMEs in the ~ 2 - 20 R© field of view us- 
ing SECCHI/STEREO data validate this assumption ( e.g., 
Poomvises et al, 2010 ). The self-similar assumption means that 
the radius of the R(t) of the flux rope is related to its heliocentric 
height H(t) by 



R(t) ^ R(T) 
H(t) H(T) 



(15) 



where H(T), the heliocentric height at time T, is = 1 AU by defi- 
nition, and R(T) is the measured radius of the magnetic cloud at 
the Earth. 

As mentioned earlier, we consider a 2-stage velocity profile 
for CME propagation, expressed by Eqs ([1} and (l2); this is sub- 
stantially different from the constant speed profile adopted in 
S09. 

Using Eqs (fTTl i. (fTSl l and (fl4l i in ( fTOb . we get the following 
expression for the total number of protons inside the CME when 
it arrives at the Earth: 



Ui 



27rL(t)R(t)Di 



Na 



^RL(t)' 



0.33 



dt. 



(16) 



The cosmic ray density inside the CME when it arrives at the 
Earth is 



A^i 



f/i 



;rR(T)2L(T) 



(17) 



where L(T) and R(T) are the length and cross-sectional radius 
of the CME respectively at time T, when it reaches the Earth. 
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Fig. 4. The muon flux along the nine directions is shown for the 
Forbush decrease on 2001 November 24. The fluxes are shown 
as percentage deviations from mean values. The solid black lines 
show the data after applying a low-pass filter (S09). The blue 
dashed line in the first panel shows the magnetic field observed 
in-situ by spacecraft. The magnetic field data are inverted (i.e., 
magnetic field peaks appear as troughs) and are scaled to fit in 
the panel. The red dotted line in middle panel shows Tibet neu- 
tron monitor data scaled down by a factor 3 to fit in the panel. 

When the CME arrives at the Earth, the relative difference be- 
tween the cosmic ray density inside the CME and the ambient 
environment is manifested as the Forbush decrease, whose mag- 
nitude M can be written as 



M 



Na-Nj 

N„ 



AN 

n7 



(18) 



We compare the value of the Forbush decrease magnitude M 
predicted by Eq ( fTSl l with observations in §3]. 

An example of an FD event observed in all 9 bins of 
GRAPES-3 is shown in Figure |4] The x-axis is the time in days 
starting from 1 November 2001 and the y-axis gives the the per- 
centage deviation of the muon flux from the pre-event mean. The 
magnitude M of the FD for a given rigidity bin is the difference 
between the pre-event cosmic ray intensity and the intensity at 
the minimum of the Forbush decrease. 
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Fig. 5. Interplanetary magnetic field and solar wind speed from 
the day 24 November 2001, The shock sheath thickness is com- 
puted by multiplying the time interval inside the dotted lines by 
the solar wind speed 



the ACE and WIND spacecrafts and estimate the shock sheath 
thickness Lshock to be the spatial extent of the magnetic field en- 
hancement. An example is shown in Figure |5] 

In computing Dj^" and D^^*'"'^'', we need to use different val- 
ues for the proton rigidity p for the ambient medium and in the 
shock sheatii; they are related to the proton rigidity Rg by 



P 



Rg 



Bo"Lsl,ock 
^shock _ Rg 



n shock J 

"0 i^slmck 



(20) 
(21) 



where BJJ is the ambient magnetic field, B^'""'* is the magnetic 
field inside the shock sheath. 



3.3. Shock-only model 

In this approach we assume that the FD is caused exclusively due 
to the shock, which is modeled as a propagating diff'usive bar- 
rier. The expression for the magnitude of the Forbush decrease 
according to this model is (Wibberenz et al 1998) 



Ms 



Ua-U 



shock 



Ua 



AU 

uT 



V™L 



shock 



Dx^ \D^ 



Dx^ 



-1 



(19) 



where Ua is the ambient cosmic ray density and Ushock is 
that inside the shock, D^^* is the ambient isotropic perpendicular 
diffusion coefficient and 0^'*^°'^'' is that inside the shock, Vs„ is 
the solar wind velocity and Lshock is the shock sheath thickness. 
For each shock event, we examine the magnetic field data from 



4. Results 

In this section we first describe various parameters needed for 
the CME-only and the shock-only models that are derived from 
observations. Using these parameters, we examine whether the 
notion of cosmic ray diffusion is valid for each model. Using 
the observationally determined parameters, we then obtain the 
(CME-only and shock-only) model that best reproduces the ob- 
served FD magnitude in each rigidity bin. 

4. 1 . Observationally derived parameters 



Table IA.2I contains the observationally determined parameters 
for each of the CMEs and their corresponding shocks in the final 
shortHst (Table[B. 
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Fig. 6. The * symbols show the Forbush decrease magnitude observed with GRAPES-3. The dashed Hne is obtained using the 
CME-only model 



The quantity "First obs" denotes the time (in UT) when the 
CME was first observed in the LASCO FOV, while RfiKt is the 
distance (in units of R©) at which CME was first observed in 
LASCO FOV and Riast is the distance at which the CME was 
last observed in LASCO FOV. The quantity Vgxp is the speed of 
CME at Riast (in km s ') and ai is the acceleration of CME in the 
LASCO FOV (in m s"^). The quantity Cd is the (constant) di- 
mensionless drag coefficient used for the velocity profile (§ 12. 2b . 
The quantities MC start and MC end denote the start and end 
times of the magnetic cloud in UT. The quantity V^^ is solar 
wind speed at the Earth (in km s ') just ahead of the arrival of 
the magnetic cloud, and Rmc is the radius of the magnetic cloud 
(in km). The quantity Bmc is the peak magnetic field inside the 
magnetic cloud (in nT). The quantity Ttotai is the Sun-Earth travel 



time (in hours) taken by the CME to travel from Sun to Earth. 
The quantity Shock arrival denotes the time (in UT) when the 
shock is detected near the Earth. The quantities B" and 3^'*°'^'' 
represent the magnetic fields (in nT) in the ambient solar wind 
and inside shock sheath region respectively (see Fig |5] for an 
example). The quantity V^^"'^'' represents the near-Eaith shock 
speed in km s"'. 

Table lA.ll contains details of the Forbush decreases associ- 
ated with each of the CMEs in the final shortlist. The magnitude 
of the FD in a given rigidity bin is the difference between the 
pre-event intensity of the cosmic rays and the intensity at the 
minimum of the Forbush decrease. It also contains the onset and 
the time of minimum, and the magnitude of the decrease in each 
bin (together with the corresponding cutoff rigidity) for each of 
the FD events in our final shortlist. 
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Table 2. Minimum ^^values for the CME-only model fits to 
GRAPES-3 data 



Event 


x' 


11 April 2001 


11.1 


17 August 2001 


1.96 


24 November 2001 


2.46 


7 September 2002 


25.5 


20 November 2003 


3.66 


26 July 2004 


27.5 



4.2. Fitting tine CME-oniy and siiocl<-oniy modeis to 
muiti -rigidity FD data 



Using the observational parameters listed in Table IA.2I we have 
computed the magnitude of the Forbush decrease using the 
CME-only (§ 13. 2b and shock-only models (§ 13. 3b . The only free 
parameter in our model is the ratio of the energy density in the 
random magnetic fields to that in the large scale magnetic field 
cr^ = (Bturb^/Bo~). Figure |6] shows the best fits of the CME-only 
model to the multi-rigidity data. The only free parameter in the 
model is cr^ = (Bturb^/Bo ) , and the best fit is chosen by min- 
imizing the X' with respect to cr~. For each FD event, the * 
symbols denote the observed FD magnitude for a given rigid- 
ity bin. The dashed line denotes the FD magnitude predicted by 
the CME-only model. We define the chi-square statistic as 



2 V (E.-Di)^ 
■^ varj 



(22) 



where Ej is the value predicted by the theoretical model Di is 
the corresponding GRAPES-3 data point and vari is the variance 
for the corresponding data points. The x^ values obtained after 
minimizing with respect to en are listed in Table |2] 

The entries in the second column ctmc in Table [3] denote the 
square roots of the turbulence parameter en that we have used 
for the model fits for each event. These values represent the 
level of turbulence in the sheath region immediately ahead of 
the CME, through which the cosmic rays must traverse in or- 
der to diffuse into the CME. By comparison, the value of cr for 
the quiescent solar wind ranges from 6-15% (Spangler 2002). 
Evidently, the CME-only model implies that the sheath region 
ahead of the CME is only a little more turbulent than the qui- 
escent solar wind, except for the 26 July 2004 event, where the 
speed of CME at the Earth was much higher than that for the 
other events. 

We have carried out a similar exercise for the shock-only 
model (§ 13. 3b . For each event, we have used the observationally 
obtained parameters pertaining to the shock listed in Table IA.2I 
Since this model needs the turbulence levels in both the ambient 
medium as well as the shock sheath region to be specified, we 
have assumed that the turbulence level inside the shock sheath 
region is twice that in the ambient medium. We find that it is not 
possible to fit the shock-only model to the multi-rigidity data us- 
ing values for the turbulence parameter that are reasonably close 
to that in the quiescent solar wind. For each event, last column la- 
beled crsj,o(.]j in Table[3]denotes the turbulence level in the shock 
sheath region that are required to obtain a reasonable fit to the 
data. Clearly, these values are an order of magnitude higher than 
those observed in the quiet solar wind. 



4.3. rL/RcME for CME-oniy modei 

Kubo & Shimazu (2010) have simulated the process of cos- 
mic ray diffusion into an ideal flux rope CME in the pres- 
ence of MHD turbulence. They find that, if the quantity 
fo(t) s RL(t)/R(t) is small, cosmic ray penetration into the flux 
rope is dominated by diffusion via turbulent irregularities. Other 
effects such as gradient drift due to the curvature of the magnetic 
field are unimportant under these conditions. Figure Qshows the 
quantity fo(t) s RL(t)/R(t) for 12 and 24 GV protons, for each 
of the CMEs in our final shortlist (Table [1}. The Larmor radius 
RlCO is defined by Eqs Q and ([T2]i and the CME radius R(t) 
is defined in Eq (fTsT l. Clearly, fo «: 1 all through the Sun-Earth 
passage of the CMEs, and this means that the role of MHD tur- 
bulence in aiding penetration of cosmic rays into the flux rope 
structure is expected to be important. 

5. Summary 

Our main aim in this work is to determine whether Forbush de- 
creases due to cosmic rays of rigidities ranging from 14 to 24 
GV are caused primarily by the CME, or by the shock associated 
with it. We examine this question in the context of multi-rigidity 
Forbush decrease data from the GRAPES-3 instrument. We use 
a carefully selected sample of FD events from GRAPES-3 that 
are associated with both CMEs and shocks. 

We consider two models: the CME-only model (§ 13.2b and 
the shock-only model (§ 13.3b . In the CME-only model, we en- 
visage the CME as an expanding bubble bounded by large-scale 
magnetic fields. The CME starts out from near the Sun with prac- 
tically no high energy cosmic rays inside it. As it travel towards 
the Earth, high energy cosmic rays diffuse into the CME across 
the large-scale magnetic fields bounding it. The diffusion coef- 
ficient is a function of the rigidity of the cosmic ray particles 
as well as the level of MHD turbulence in the vicinity of the 
CME (the sheath region). Despite the progressive diffusion of 
cosmic rays into it, the cosmic ray density inside the CME is 
still lower than the ambient density when it reaches the Earth. 
When the CME engulfs the Earth, this density difference causes 
the Forbush decrease observed by cosmic ray detectors. In the 
shock-only model, we consider the shock as a propagating dif- 
fusive barrier. It acts as an umbrella against cosmic rays, and 
the cosmic ray density behind the "umbrella" is lower than that 
ahead of it. 

We have obtained a list of Forbush decrease events observed 
by the GRAPES-3 instrument using the shortlisting criteria de- 
scribed in § |2] For each of these shortlisted events, we have used 
observationally derived parameters listed in Table [A2l for both 
the models. The only free parameter was the level of MHD tur- 
bulence (defined as the square root of the energy density in the 
turbulent magnetic fluctuations to that in the large-scale mag- 
netic field) in the sheath region. 

6. Conclusions 

Figure |6] shows the results of the CME-only model fits to multi- 
rigidity data for each of the shortlisted events. For the shock-only 
model, we use the turbulence level in the shock sheath region as 
the free parameter. Table |3] summarizes the values of these tur- 
bulence levels that we have used for each of the FD events in 
the final shortlist. These values may be compared with the es- 
timate of 6-15 % for the turbulence level in the quiescent so- 
lar wind (Spangler 2002). We thus find that a good model fit 
using the CME-only model requires a turbulence level in the 
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Fig. 7. The quantity /o s Rl/R versus distance as a CME propagates from the Sun to the Earth. The dashed and continuous lines 
represent 24 GV and 12 GV protons, respectively. 



Table 3. Turbulence levels in the sheath region required by the 
models 



Event 


CME-only model 


shock-only model 




0"MC 


O'Shock 


11 April 2001 


9.4% 


100% 


17 August 2001 


13% 


180% 


24 November 2001 


28% 


400% 


7 September 2002 


13% 


100% 


20 November 2003 


6.7% 


400% 


26 July 2004 


46% 


200% 



sheath region that is typically only a little higher than that in 
the quiet solar wind, which is generally consistent with observa- 
tions. On the other hand, a good fit using the shock-only model 
demands a turbulence level in the shock sheath region that is 
often an order of magnitude higher than that in the quiet solar 
wind, which is somewhat unrealistic. The results summarized in 
Table |3]imply that, for FDs involving protons of rigidities rang- 
ing from 14 to 24 GV, the CME-only model is a viable one, 
while the shock-only model is not. Given the remarkably good 
fits to multi-rigidity data (Figure |6] the reasonable turbulence 
levels in the sheath region demanded by the CME-only model 
(Table [3) and because the FD minima usually occur well within 
the magnetic cloud (Table IaTTT i. we conclude that CMEs are the 
dominant contributors to the FDs observed by the GRAPES-3 
experiment. 
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Table A.l. Derived parameters of Forbush decrease for events, 



11 April 2001 

' FD onset times are in UT, 1 1 April 2001 
^ FD minimum times are in UT, 12 April 2001 



Shock 


Arrival 


: 11 April 2001, 14:0f 














Magnetic Cloud 


Start 


: 11 April 2001, 23:0C 


) 




End 


: 12 April 2001, 


18:00 


Quantity 


NW 


N 


NE 


W 


V 


E 


SW 


S 


SE 


FD Magnitude (%) 


1.02 


1.12 


0.91 


1.38 


1.40 


1.05 


1.36 


1.29 


0.93 


Cut-off Rigidity (GK) 


15.5 


18.7 


24.0 


14.3 


17.2 


22.4 


14.4 


17.6 


22.4 


FD Onset (UT) 


12:57 


13:55 


16:04 


11:16 


12:00 


12:57 


08:24 


10:20 


12:28 


FD minimum {UT) 


19:00 


19:00 


19:00 


17:00 


18:00 


18:00 


14:00 


15:00 


16:00 



17 August 2001 

' FD onset times are in UT, 17 August 2001 

^ The * FD onset times are in UT, 16 August 2001 

^ FD minimum times are in UT, 18 August 2001 



Shock 


Arrival 


: 17 August 2001, 


11:00 










Magnetic Cloud 


Start 


: 18 Aug 


ust2001. 


00:00 




End 


: 18 August 2001 


21:30 


Quantity 


NW 


N 


NE 


W 


V 


E 


SW S 


SE 


FD Magnitude (%) 


1.17 


1.12 


0.83 


0.97 


1.03 


0.79 


1.12 1.07 


0.80 


Cut-off Rigidity (GK) 


15.5 


18.7 


24.0 


14.3 


17.2 


22.4 


14.4 17.6 


22.4 


FD Onset (UT) 


04:19 


01:55 


01:12 


23:17* 


22:34* 


22:05* 


00:00 23:31 * 


23:17* 


FD minimum {UT) 


4:00 


3:00 


2:00 


6:00 


5:00 


5:00 


23:00 22:00 


20:00 



24 November 2001 

' FD onset times are in UT, 24 November 2001 
^ FD minimum times are in UT, 25 November 2001 



Shock 


Arrival 


: 24 November 2001 


06:00 












Magnetic Cloud 


Start 


: 24 November 2001 


17:00 




End 


: 25 November 2001 


, 13:00 


Quantity 


NW 


N 


NE 


W 


V 


E 


SW 


S 


SE 


FD Magnitude (%) 


1.42 


1.34 


0.94 


1.67 


1.56 


1.16 


1.34 


1.36 


1.10 


Cut-off Rigidity (GV) 


15.5 


18.7 


24.0 


14.3 


17.2 


22.4 


14.4 


17.6 


22.4 


FD Onset {UT) 


03:07 


03:21 


03:07 


04:05 


03:21 


02:52 


04:05 


03:07 


01:24 


FD minimum {UT) 


17:00 


15:00 


16:00 


14:00 


15:00 


17:00 


15:00 


16:00 


19:00 



7 September 2002 

' FD onset times are in UT, 7 September 2002 
^ FD minimum times are in UT, 8 September 2002 



Shock 


Arrival 


: 7 September 2002, 


14:20 












Magnetic Cloud 


Start 


: 7 September 2002, 


17:00 




End 


: 8 September 2002 


16:30 


Quantity 


NW 


N 


NE 


W 


V 


E 


SW 


S 


SE 


FD Magnitude (%) 


0.41 


0.59 


0.63 


0.89 


0.97 


0.83 


1.08 


1.07 


0.84 


Cut-off Rigidity (GK) 


15.5 


18.7 


24.0 


14.3 


17.2 


22.4 


14.4 


17.6 


22.4 


FD Onset {UT) 


17:03 


17:17 


16:05 


14:52 


14:52 


15:07 


15:50 


16:19 


17:03 


FD minimum {UT) 


15:00 


14:00 


14:00 


13:00 


13:00 


14:00 


15:00 


16:00 


16:00 



20 November 2003 

' FD onset times are in UT, 20 November 2003 
^ The * FD onset times are in UT, 19 November 2003 
^ The ** FD onset times are in UT, 21 November 2003 
'* FD minimum times are in UT, 24 November 2003 



Shock 


Arrival 


: 20 November 2003, 


07:30 












Magnetic Cloud 


Start 


: 20 November 2003, 


10:06 




End 


: 21 November 2003 


, 00:24 


Quantity 


NW 


N 


NE 


W 


V 


E 


SW 


S 


SE 


FD Magnitude (%) 


0.95 


0.93 


0.81 


1.19 


1.16 


0.97 


1.17 


1.20 


0.93 


Cut-off Rigidity (GV) 


15.5 


18.7 


24.0 


14.3 


17.2 


22.4 


14.4 


17.6 


22.4 


FD Onset {UT) 


21:22* 


08:10** 


02:53** 


01:55 


10:48 


20:38 


06:58 


10:19 


15:07 


FD minimum {UT) 


5:00 


4:00 


2:00 


4:00 


4:00 


4:00 


3:00 


3:00 


2:00 



26 July 2004 

' FD onset times are in UT, 26 July 2004 
^ FD minimum times are in UT, 27 July 2004 



Shock 


Arrival 


: 26 July 2004, 22:20 












Magnetic Cloud 


Start 


: 27 July 2004, 02:00 






End 


: 27 July 2004, 24:00 


Quantity 


NW 


N 


NE 


W 


V 


E 


SW 


S SE 


FD Magnitude (%) 


2.06 


2.04 


1.67 


2.40 


2.12 


1.62. 


1.77 


1.67 1.40 


Cut-off Rigidity (GV) 


15.5 


18.7 


24.0 


14.3 


17.2 


22.4 


14.4 


17.6 22.4 


FD Onset {UT) 


14:24 


15:22 


17:17 


14:10 


15:36 


18:00 


16:05 


18:29 20:38 


FD minimum {UT) 


10:00 


11:00 


13:00 


10:00 


11:00 


13:00 


12:00 


14:00 15:00 



Notes. For each event, first row: magnitude of Forbush decrease (FD magnitude), second row: cut-off rigidity, third row: Forbush decrease onset 
time, fourth row: FD minimum time. The units in which each quantity is expressed is given in parentheses in the first column. 
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Arun Babu et al.: High rigidity Forbush decreases 
Table A.2. Observed parameters of CME & Shock for different events 



Event 


11 Apr 2001 


17 Aug 2001 


24 Nov 2001 


7 Sep 2002 


20 Nov 2003 


26 Jul 2004 


CME details 


" First obs. ( UT) 


10/04/01,05:30 


15/08/01, 23:54 


22/11/01,22:48 


5/09/02, 16:54 


18/11/03,8:50 


25/07/04, 14:54 


Rflrst (^o) 


2.84 


3.38 


4.77 


4.12 


6.3 


4.22 


'Rlas. (Re) 


18.1 


25.9 


25.9 


17.0 


27.5 


21.9 


"VexpC/cm.-') 


2880 


1410 


1370 


1860 


1650 


1370 


'' a, {m s'^-) 


211 


-31.7 


-12.8 


43.0 


-3.29 


7.0 


■'Cd 


0.325 


0.163 


0.09 


0.312 


0.333 


0.016 


s MC start (UT) 


11/04/01 ,23:00 


18/08/01,00:00 


24/11/01, 17:00 


7/09/02, 17:00 


20/11/03, 10:06 


27/07/04, 2:00 


" MC end (UT) 


12/04/01, 18:00 


18/08/01,21:30 


25/11/01, 13:00 


8/09/02, 16:30 


21/11/03,00:24 


27/07/04, 24:00 


• YZ (km s-') 


725 


600 


730 


544 


750 


900 


'Rmc (km) 


2.48 X 10' 


2.32 X 10' 


2.63 X 10' 


2.30 X 10' 


1.89 X 10' 


3.56 X 10' 


'Bmc (nT) 


34.5 


25.6 


20 


22.9 


50 


25.3 


' T.otai (hours) 


42.1 


44.9 


41.1 


48.9 


42.1 


35.7 


Siiock 


Siiock arrival (UT) 


11/04/01, 14:06 


17/08/01, 11:00 


24/11/01,6:00 


7/09/02, 14:20 


20/11/03,7:30 


26/7/04, 22:20 


"' B" {iiT) 


4.5 


5 


5 


5.8 


7 


5 


'• B^"""^" {nT) 


32.5 


33 


41.5 


23 




26.1 


"Yf^'^'ikms-') 


670 


501 


948 


550 




893 



Notes. *°' First obs. : time in UT when the CME was first observed 

'*' Rfirst : distance in (Re), where CME was first observed. 

•'■' Riast : distance in (Rq), where the CME was last observed in the LASCO FOV. 

*'" Vexp : speed of CME in (km i-"') at Ri^st 

<'■> ai : acceleration of CME in LASCO FOV in units of (in s'^) 

'-'^' Cd : the constant dimensionless drag coeificient used for the velocity profile (§ 2.2). 

'*' MC start : Magnetic cloud start time in UT. 

'•''^ MC end : Magnetic cloud end time in UT. 

''' ^sw" • solar wind speed in (km .v"') at Earth during the arrival of the magnetic cloud 

'■'' Rmc : radius of magnetic cloud in (km ) 

'*' Bmc '■ magnetic field inside the magnetic cloud (iiT) 

*'' Ttoiai : The Sun-Earth travel time in hours for the CME. 

(m) ga . juagnetic field in the ambient solar wind (in nT ) 

(n) gshock . jjjagi^gtJQ ggijj inside the shock sheath region (in nT) 

'°* V*"'^'' : near-Earth shock speed (km .?"') 



12 



